Hall conductivity in the presence of spin-orbit interaction and disorder 
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Starting from the Kubo formula, we expand the HaU conductivity using a cumulant approach 
which converges quickly at high temperatures {ksT > energy differences of initial and final scattering 
states) and can be extended to low temperatures. The theory can deal with the sign, the ordinary and 
the anomalous contributions to the Hall effect. When applied to include the spin-orbit interaction 
to first order, we recover what is essentially the Karplus-Luttinger result for the anomalous Hall 
effect. Contact is made to the Chazalviel and Nozieres-Lewiner formulae. A side-jump type formula 
is obtained by using an exact application of linear response. We show that there exists an exact rigid 
Hall current which is not a Fermi level property. We introduce a relationship between mass and 
diffusivity which allows us to generalize the theory to strong disorder and even introduce a mobility 
edge. The formalism provides a systematic and practical way of analyzing both ordinary and 
anomalous contributions to the Hall conduction including the changes of sign, and in the presence 
of serious disorder. As a byproduct of the method, we show that the anomalous Hall coefficient can 
vary with resistance to the power n, with 1 < n < 2 depending on the degree of coherence. 

PACS numbers: 72.10.Bg, 72.15.Gd, 75.47.-m, 72.80.Ng, 71.70.Ej 
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The Hall conductivity of materials exhibits a wide 
and rich variety of behavior. The interpretation is, in 
general, still very difficult, even though, in principle, 
the information is contained in the Kubo formula. 
This is true at least in linear response to the applied 
electric field. The Kubo formula is therefore the starting 
point of our analysis of Hall conductivity. We include 
the spin-orbit interaction and discuss the so called 
anomalous Hall effect (AHE). Previously, most work 
was focused on understanding the sign change of the 
Hall coefficient in ordered and disordered systems^, 
localization, and the Quantum Hall effect (QHE)^. 
More recently, the problem has been to understand the 
effect of magnetism and of many body corrections on 
the Hall effect. Considerable progress has been made 
recently by integrating the spin-orbit coupling into the 
Bloch wavefunction formalism of Kane in crystals, and 
applying these wavefunctions and other first principle 
numerical methods^ to study the AHE effect in magnetic 
materials^. In these papers, the emphasis is on order, 
and disorder is only represented by a uniform lifetime. 
Most magnetically doped semiconductors and magnetic 
alloys can, however, not really be considered to be 
in the weak scattering regime, and therefore, in the 
present approach, we have inverted this priority. We 
emphasize the absence of Bloch symmetry rather than 
its presence 
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The aim is to develop practical formulae which 
can deal with disorder, the sign, the side jump Hall 
effect problem, and the Skew scattering/intrinsic Hall 



The AHE is now a well established observation in 
magnets, and a number of explanations have been 
proposed^, and will be discussed here. All explanations 
are based, in one way or another, on the spin-orbit 
mechanism. Thus, it is accepted that spin-orbit coupling 
causes the anomalous Hall conductivity contribution. 

Much progress has been made recently on under- 
standing the origin of the AHE and this has generated 
exciting new physics. Traditionally it was thought that 
only the spin-orbit Skew scattering mechanism^'- i^^^^'^ , 
gives rise to a magnetism (M^) dependent Hall co- 
efficient. This process depends on the conductivity 
relaxation time or relaxation time squared, depending 
on whether the skew scattering is itself rate determining 
for conductivity or not. Traditionally, skew scattering is 
derived as an extrinsic effect, i.e. it is due to impurity 
scattering, and not to the host spin-orbit interaction. 
The intrinsic spin-orbit interaction produced by the 
host crystal potential had been invoked as a source of 
AHE by Karplus and Luttinger^ but later rejected by 
Smit^^ who claimed that the intrinsic effect is negligible. 
Earlier, it had been shown by Mott and Masseyi^ 
that electron scattering from Coulomb potentials is 
asymmetric with respect to spin direction, with spin up 
going more to one side, and those with spin down more 
to the other. When the electron gas is magnetized, there 
is a net transverse Hall current. When the asymmetric 
scattering is at impurity sites, this process is called skew 
scattering as mentioned above, and has been discussed 
by several authorsiiii^. The important point about 
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skew scattering is that the net spin-orbit coupling at 
the impurity site, can be strongly enhanced by the 
host crystal. This is really what makes this process so 
important. For a detailed account of the history and 
progress in understanding the AHE see Ref. i8|. 



II. THE KUBO FORMULA FOR THE 
LONGITUDINAL AND TRANSVERSE 
CONDUCTIVITY 

The frequency (w) dependant conductivity in linear 
response to an electric field is usually written asF^^t^ 



ihe -r-^ {( 



\v^\P){(3\v,\a) /(£o)-/(£/3) 

— 6(3 + Hu + iS Sjj — Sa 



(1) 

This form is general for any exact set of eigenstates 
\a) and energies Sa, spin summation is implied. The 
/(e) are the Fermi functions and SI is the volume. The 
velocity operators are given by Heisenberg's equation 
of motion. 

For the general case of a material which need not 
be periodic, the spin orbit term in the Hamiltonian is 



Hfin — 



(VV{v) X p) . ^, 



where 



V{v) = J2 



eZn 



ineeolr - R„| 



eFx. 



(2) 



(3) 



In Eqs. ^ and ([3]), m is the bare electron mass, c is 
the speed of light, p is the momentum operator and cr 
is the Pauli spin operator which is a vector containing 
the Pauh's matrices i.e. [(Jx,cry,az]- In Eq. ([3]), F is 
the external applied electric field, e the electric charge, 
Zn is the effective local charge, eeo is the permittivity, r 
and X are positions operators for the charges and R„ is 
the position of the fixed ions that make the lattice. The 
velocity operators can then be written as 



(4) 



y 



eB,x 



V,V{r)a,~V,V{r)ax , (5) 



where we have used the Landau gauge for the vector 
potential, A = {0,B,x,0). The v° and v° are 

respectively. The choice of the minus sign 



and 



ih d 
m dy 



in —eBzx/m implies that e = — |e| since the kinetic 
part of the Hamiltonian for a charge g in a field is 
T = l/2m(p-gA)2. 

The spin dependent terms in Eqs. @ and ^ can 
be important in magnets. We shall consider them 



explicitly in Section IVl 

The Hall effect is given by the antisymmetric part 
of the transverse conductivity'iiiii. To write the 
antisymmetric Hall conductivity, we note that for 
non-zero frequency it is the anti-hermitian partii 
a'^y = ^{(Jxy — cTyx)- the dc limit we have 



he' 



lim ^ -i{a\vx\f3){P\vy\a) 



a, 13 



{ef3 — Ea) 



S' 



(6) 

We shall now derive the Hall conductivity in a general 
way which will include both the skew scattering and the 
intrinsic contributions. 



III. THE HALL EFFECT USING A CUMULANT 
EXPANSION OF THE KUBO FORMULA 

One way to derive the contributions to the Hall effect, 
with disorder present, is to rewrite the Kubo formula 
using the Heisenberg equation of motion: 



{a\x\l3) = —ih 



£a- ' 



(7) 



which is true in any finite box (length L) without dissi- 
pation. Then Eq. ([6]) becomes 



a. 



xy 

— lim > 1 



a,/3 



:\P){/3\vy 



(8) 



We drop the a superscript for the Hall conductivity, but 
it will be implicit that we refer to the antisymmetric part 
unless otherwise mentioned. In order to demonstrate the 
cumulant technique, we first consider the limits of high 
temperatures when {sf3 — < ksT and, of small mag- 
netic fields kgT ^ hujc- This is the limit when the ma- 
trix elements are dominated by intraband scattering with 
weak B-field. We expand the function 



£0 — £a 



e J — e s 



knT 



t- ^ ( £f3 — £a 



2! V knT 



e ''B^ 



(9) 



and substitute Eq. jO]) and Eq. ^ into Eq. (P. This 
gives us a cumulant expansion of the Hall conductivity 
in powers of -jp^ where one can use the operator identity 



{£f3-£a)Aaf3 = {a\[A,H]\[3), 



(10) 



where [A, H] is the commutator of A with the Hamilto- 
nian H and A^p = {a\A\(3), to reduce and evaluate the 
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terms generated by the expansion. 

The high T expansion converges rapidly as soon as 
fcsT is larger than the typical energy difference in the 
matrix elements. In particular, this is obviously true 
as soon as ksT exceeds the bandwidth. But, we can 
also have a more general expansion that may also be 
valid at low temperature though the justification in this 
case is more complicated. To obtain the more general 
cumulant expansion at any temperature, we observe 
that in Eq. ^ the dominant region is around ~ Eq. 
Thus, we carry out a complete Taylor expansion of /(s/j) 
around Sa and obtain 



/(ea) - f{ep) 



(11) 



When we substitute Eq. (fTT|) (stricly equivalent to Eq. ([9|) 
at high T) into the Kubo formula (Eq. ([8])), we obtain 
the first order term 



r{l} 
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a,/3 



df{ea)\ {sfj 



dSa J {Sfj - £a)2 + 

x{a\x\P) {(3\vy\a) 



lim > 

a. 13 



dSa 



{a\xvy \a) 
5^ 



{a\x\p){l3\vy 



and therefore are lower order corrections. From the 
mathematical structure and to keep the terms linear in 
Bz, it is thus essential therefore only to keep the second 
order cumulant at low temperatures. This will now be 
shown step by step as we proceed with the analysis of 
the various contributions to the Hall current. 

Let us also remember that in Eq. HH), \a) is a 
magnetic field dependent exact eigenstate. 

In order to rewrite the matrix element in a simple 
way, we show in Appendix \^ that we can write the 
derivative of the eigenvalue with respect to the magnetic 
field as 



dBz 



-e {a \xvy \ a) 



a 



a), (13) 



where Vy does not contains the spin-orbit term (Eq. ([5])) 



contrary to Vy in Eq. (|12p . The spin-orbit term of Vy 
will be treated later, in Section IV CI 

The Hall conductivity only comes from the first 
term on the RHS of Eq. the orbital term. The 

way to handle this is to introduce initially two different 
magnetic fields, one acting on the orbital part Borb and 
one giving the Zeeman energies B^ ■ We can thus rewrite 
Eq. nil) as 



= 4i:(-^)<°|"'»i")- m 

where we have used the identity lim^^o „g )'^+s'^ ~ 
TTS{ea — sp) lima^o S — 0. 

Apart from this first order term (Eq. (fT2|) ). we ob- 
tain an infinite series of higher order cumulant which 
can be generated via Eq. (fTO|) . These higher order 
cumulants represent correction to the first order result. 
The first order term will turn out to have very simple 
physical interpretation. The corrections generated by 
the higher order cumulants resulting from Eq. ([ni), and 
when evaluated at low temperatures, must be studied 
in detail as it is not immediately self-evident that they 
represent lower order corrections. This is done later in 
the paper, but at this stage, it is already possible to note 
that the corrections represent, via Heisenberg's equation 
of motion, and from Eq. (jlOp . higher and higher time 
derivatives of the velocity operator. When using the 
effective mass Hamiltonian, one can show that beyond 
the second order cumulant, the corrections which are 
linear in the cyclotron frequency or magnetisation all 
scale with the disorder potential and spin-orbit coupling, 



ne^\ dec. J 



Only the derivative of the energy eigenvalues with re- 
spect to the orbital field Borb gives the Hall conductivity. 
An advantage of this representation is that the spin-orbit 
energy can now be treated in first order perturbation 
theory, as we shall see later. 

The way disorder should be treated for numerical 
calculations is as follows. The Hall conductivity as 
given by the general conductivity formula (Eq. ([T])) or 
the formulae obtained using the cumulant approach 
(Eq. (|14p and equations thereafter) is evaluated for a 
specific configuration of disorder. The wavefunctions \a) 
and energies are therefore the exact corresponding 
eigenstates and eigenvalues. Now the procedure is 
repeated for all the possible configurations and averaged 
with the appropriate weighting factors. 

Before examining the higher order terms (Appen- 
dices |B] and |C]) , let us understand the significance of this 
result and compare with other well known approaches. 
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IV. COMPARISON TO OTHER THEORIES 



where N 



p{e f)'m*v'j. 



A. The Streda result 



We can rewrite Eq. p4|) in the form 



r{l} 



(15) 



E=ef 



where p{E',Bz) is the density of states with magnetic 
field. 

This expression is the quantum term of the Streda 
formulaii^. This author reduced the Kubo formula to 
two terms called ct^^ and crly (Eq. (fTSjl above). The 
antisymmetric part of the other term, <J^y, should be 
contained in the remainder of the cumulant expansion. 
Our term (Eq. (jlSp) differs from Streda'si^ by a minus 
sign. The sign problem can be traced in Ref. ^^ to one 
transformation (Eq.(ll) of Ref. ITsl ) where there should 
be a minus sign on the RHS. 



V. THE INFLUENCE OF THE SPIN ORBIT 
COUPLING 

Spin-orbit coupling introduces a number of new con- 
tributions to the Hall current. New terms arise, due to 
the spin dependent velocities from Eqs. ([U and ([5]), and 
from the effect of the spin-orbit interaction on the energy 
levels. We examine this last effect first. Also, as pointed 
out in Section IIVB( for the class of problems where we 
have a periodic system -I- impurities, wherever the bare 
mass appears in the following text (except when it comes 
from the spin-orbit Hamiltonian) , we can replace bare 
mass with the effective mass (m*) and drop the periodic 
part of the Hamiltonian in the remaining analysis. 



A. The effect of the spin orbit interaction on the 
Hall current from the changes in the energy levels 
and wavefunctions 



B. The classical limit 



If, instead of Eq. (O, we write as in Ref. [Ti 



{a\x\p) 



-ih- 



{a\v.\P} 

^ — ^13 + - 



(16) 



where r is a lifetime, we can reduce the first term of the 
cumulant expansion to the classical result. By adding 
a lifetime, we effectively assume that the electrons are 
subject to resistive scattering processes. Using the 
effective mass Hamiltonian for the periodic part of 
the Hamiltonian also allows us to replace m with the 
effective m*. The finite lifetime then represents all 
the scattering processes that break the translational 
invariance of a Bloch electron with effective mass to*. 
This includes disorder and electron-phonon scattering 
treated in the Born approximation. When working in 
this approximation, there is no configuration average to 
be done anymore. 

We consider only the term linear in Bz in Eq. ([8|). 
the one involving the eBzx/m part of Vy (see Eq. ([5])), 
which gives the diagonal mass tensor term. Then, 
Eq. ([5]) together with Eq. (fTB|) yields the cumulant. 



2 ^ 



{a\x'^\a) 



dfi^a) \ sr-^ 2 1/ I I m |2 



(17) 



This leads to the very well known classical result 
{1} A^eV eBzT eB^T 



Consider the first order cumulant result with the zero 
order velocity operators (including the —^^^ term). The 
advantage of the cumulant expansion is that it allows 
us to analyze a very complex phenomenon, the effect of 
the spin-orbit coupling on the Hall conductivity, via the 
eigenstates. We write, to first order in perturbation the- 
ory, 



ec.{Bz)^el{Bz) + {a\Hso\a) 



(19) 



The spin orbit Hamiltonian being dependent upon the 
Pauli's matrices, we should remember that the state \a) 
must now be a spinor (two components vector). The 
action of taking the bracket will leave a scalar and, as we 
use the same state at this order, only the z component 
survive. 



dea{Bz 

dBorb 



= —e{a\xv,,\a} 



del{Bz 

dBorb 



where 



An — 



dBo 



heZ„ 



a 



Am?c^{A:TTe£Q) |r — R„ 



(20) 



(21) 



a" = (ajcTzIa) and li^z is the z component of the orbital 
angular momentum at site i. 

The first term in Eq. (|20p is the one we examined 
above in the classical limit and is intuitively very attrac- 
tive. The magnitude of the Hall current per eigenstate 
is related to the sensitivity of the energy level to an 
external magnetic field. Its sign depends on whether 
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the magnetic field increases or decreases the energy of 
the eigenstate. In particular, it also follows that the 
contribution of a localized state is negligible. In reality, 
localized states should actually give exactly zero. It 
trivially follows that the zero is recovered only after 
summing the remaining contributions in the cumulant 
series. Keeping only the first cumulant does not give 
the exact result when the level in question is a localized 
level, with discrete energy levels. But the first order 
result is close to zero, and therefore can be said to 
represent a good approximation. See Appendix [b1 and ICl 
for the analysis of the higher order cumulants. 

For dclocalized states, there is more information in 
the first term of Eq. (j20|l . Normally, for weak scattering, 
when e f is near the top of the band, we have the hole 
sign, because the magnetic field can only lower the 
energy near the top of the band. Near the bottom of a 
band, we have the electron sign because the magnetic 
field confines the carrier and raises the energy of the 
electrons. This rule is also true for disordered eigen- 
states. The effect of the magnetic field on the energy 
levels can be evaluated in second order perturbation 
theory in the presence of disorder . 

The anomalous term can now be studied by going 
to first order perturbation theory in _B-ficld with exact 
eigenstates. The term — /i^L • Borb generates 

H = |ao) + Y.^l3o\ - A^sL • Borb\ao) ■ (22) 

Pa 

Substituting Eq. ^ in the second term of Eq. pl)l . 
keeping only the z-components, we obtain terms in the 
energy which involve a factor of the type 

^ jJ,Z 

Ag^^^E (23) 

The sign of the anomalous process depends upon the sign 
of a quantity which is closely related to the electron g- 
shift and which itself can be electron-like or hole-like. 
Thus the first order cumulant can be written as 

Let us compare the relative magnitude of the two terms 
of Eq. dHl). At very low T ( -|| = S{e - Sf) ) we can 
rewrite the Hall conductivity (Eq. and the current 

is 

Jy^4l^^"^P^+e'p{ef)^(a.),,Af^ef)F,. (25) 

With p{ef) ~ 10^5 /m^J and = 10* V/m we have for 
the anomalous contribution J°" = 10'^Ag{az)Fx A/m^, 
which is the same order of magnitude as the normal 
process with TV = 10^^/m'^, giving the normal Hall 



current = 10'' A/m^ or 10^{huJcT)Fx A/m^ where 
hwcT 10""^. In principle, the normal and anomalous 
terms can have opposite signs. The results of Eq. ([24|) 
and Eq. (|25p are very elegant, and by writing Ag = 2 — g* 
we have essentially recovered the Fermi level version 
of the Chazalviel22 and Nozieres and Lewiner— result. 
Chazalviel"-- computes the single carrier wave packet 
motion in an electric field without using the Kubo linear 
response formalism. We also note that in this form, 
the anomalous term, apparently has no dependence on 
the relaxation time. Finally, and most importantly, 
it can also be related to the Karplus and Luttinger— 
result in the limit of high temperature when the Karplus 
and Luttinger energy gap between the Kane Luttinger 
subbands is taken as A ~ kgT and Ag ~ 

In the Karplus-Lutttinger— Bloch wavefunction 
formalism, the anomalous Hall effect, even though it is 
intrinsic, is a Fermi level property at low temperatures. 
In Rcf. 10, the spin-orbit interaction is treated in first 
order perturbation theory using Bloch functions. We 
have also used first order perturbation, but the Fermi 
level property, here, is a result of keeping the first 
order cumulant. An interesting point is that Karplus 
and Luttinger did not use the Kubo formula. They 
arrived at a similar expression, except that the matrix 
elements are always interband matrix elements. The 
reason is that their starting point is the Bloch function, 
so that in the absence of an explicit treatment of 
disorder scattering, only interband matrix elements 
are left when the matrix elements of position and 
momentum operators are considered. Note that one 
way of calculating the g-shift, when we have weak dis- 
order, is to use the Kohn-Luttinger wavefunctions. One 
can compute the f;-shift At/k.a in the exact band states^^. 

In summary, in this section, we have derived a re- 
sult which can be related to the Karplus-Luttinger 
intrinsic AHEiiS and we have made contact with the 
Chazalviel^S, Nozieres and Lewiner^-'^ and Sinovai^ 
results using a simple and unified formalism. 

B. The Side Jump Hall current 

The new velocity terms generated by the spin-orbit 
coupling in Eqs. ^ and ([5]) have been shown by Wolfle 
and Muttalib^^ to give rise to a term of a form called 
the "side jump Hall effect" in the linear response Kubo 
formula, which is 

J, = -7Ve2(a,)^F.. (26) 

However, this result is obtained only after a complex 
diagrammatic sum of potential scattering events. A 
similar result was derived previously and more simply 
from the same velocity term by Lyo and Holstein'^^ using 
scattering theory and originally by Berger^^. We shall 
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now rederive a contribution to the Hall current which 
has the same structure but is more closely related to 
the Rashba effect--. It is derived using linear response 
with the external field induced spin-orbit term in the 
Hamiltonian, below. The linear response analysis will 
show us how the presence of the lattice goes to modify 
the Hall conductivity even when we have disorder. 
Treated in linear response, the Rashba coupling term 
will give a side jump like Hall current. We do not use 
the cumulant approach or effective mass approximation 
because it is easy to treat this term exactly. 

Consider again that part of the spin-orbit coupling 
which is itself dependent on the applied external field 
(see Eqs. ([3|)). This term is a contribution to the total 
energy which directly depends on the external field. As 
part of the Hamiltonian, this term creates a departure 
from equilibrium, and must therefore be treated on the 
same footing as the usual electric potential eF^x. We 
therefore start from first principles, with the density 
matrix. When we take all such terms in the Hamiltonian 
as the perturbation i/pertj the change in the density 
matrix is given by 



AD 



Q/3 



{a\H, 



pert I 



g/3 - Eq 



(27) 



where 



H, 



pert 



eF,x + (VK.t(r) x p) • cr, (28) 



with Vext{ic) = -eFxX. 

The second term of Eq. pS)) involves the external 
applied electric field, and one can use linear response 
and ask: what Hall current does it produce in the 
presence of disorder? We can evaluate the thermally 
averaged velocities in the usual way. We consider the 
external-field-independent eigenstates, include disorder 
and the Zeeman splitting. These states can therefore be 
picked to have either spin up or spin down eigenstates 
in a chosen direction. Thus the y-current produced by 
the second term from Eq. ^7} and Eq. is given 

by (spin projection in z-direction and keeping only the 
y-momentum term) 



-F,Y.(f3\vy\a) 

a, 13 



£0 — £a 



(29) 



If we now use mvy = py and the sum rule (to be discussed 
in detail in Section IVDp 



2ma '-^ ep - Sa 



(30) 



where ji = x ov y, applied to the y-operator, we have the 
very simple and elegant result 



which apparently only depends upon the effective mass. 
This is true as long as Eq. ([50]) can be used to define 
effective mass, i.e. if all incoherence is neglected. The 
interpretation of this sum rule in the Kubo formula 
context is not trivial. If we include the entire infinite 
spectrum in the sum of Eq. (j30p , then we have the trivial 
result fUa = m, and we obtain a contribution which 
looks exactly like the Berger— and Lyo and Holstein— 
side jump Hall current but is clearly based on a different 
logic (Rashba term). Note that in our theory, the spin 
polarization is now summed over the entire band and is 
not just the fermi level spin. 

In the framework of a finite band model, one can 
interpret Eq. pop as the effective mass. The spin-orbit 
terms associated with the lattice can be included in 
Eq. (|30p . In the weak disorder limit, one could compute 
Eq. pO|) using the Kane method. One can see that, in 
the Kane model, Eq. ([50)) is indeed the effective mass. 
The effective mass correction in Eq. (j3ip can, then in 
principle, increase the current up to two orders of mag- 
nitude (InSb for example). But the side jump (Eq. ([31]) ') 
is, even with effective mass, for extended states, 
much smaller than the Karplus-Luttinger contribution. 
We have Jy — 10'"^^ {az)NFx A/m in 2-dimensions. 
In 3-dimensions, with the same numbers we have 
10^^ Fx:(az){m/m*) A/m^ compared to 10^Fx{az)Ag 
A/m^ from Eq. ([^5|) . Note that the huge Bloch function 
enhancement evaluated in various forms by ChazalvieP— , 
Fivaz?S and Berger— and which make the side jump 
term important does not appear in Eq. ([?T|) . Chazalviel, 
for example, used the Heisenberg commutator for Vy 
and Kane wavefunctions to derive the Hall velocity and 
then derives a similar expression for the Hall current 
using Drude theory. We have derived an expression 
which is similar in structure to what is called the side 
jump Hall current in the literature. The derivation we 
have used is however not the same as that of Berger—, 
Lyo and Holstein^iand Wolfie and Mutallib— . In our 
formula, the spin-orbit coupling can be enhanced by the 
lattice, but the effect is relatively small, and directly 
related to the effective mass lowering (see Eg . (I5T]) ) . The 
so-called side jump theories^li^ in which the mechanism 
is due to the spin-orbit scattering induced sideways 
jump at impurity potentials, and the associated large 
enhancements caused by Bloch functions, have not been 
recovered using the present Kubo formula method. 



To complete the analysis of Eq. ([31]) we need a 
discussion of the sum rule of Eq. ([50)) and we will defer 
this to Section |VD] because a similar term is encountered 
in the next section. 



7 



C. The effect of the spin dependent velocity on the 
Hall current: The terms which are due to the 
internal potentials 

Let us consider the contributions to the Hall current 
which results from including the contribution of the re- 
maining spin-orbit velocity terms (Eq. ([¥]) and (l5|)) in the 
Kubo formula. These now involve the lattice potentials 
as sources of velocity. For Coulomb potentials, we have 



{a\v,\P){p\vy\a) 



eZ„h 



y 



^ Am'^c'^ (4:Treeo) |r 
1 



f3){(3\vy\a) 




v){v\y'Jz\P) {f3\vy\a) 



(32) 



a) {a \ycrz\P) {P\vy\a) 



where ^ = 4^.gf("4^,,^) ■ 



Two approximations were made here. takes al- 

ternatively positive and negative values and the term 
involving it would be zero if r were not present. But, 
even if we include r, it will always give a smaller 
contribution compared to the first and we therefore 
neglect it. Second, |j is local and therefore cannot 

couple different sites but, in tight-binding for example, 
could couple different orbitals at the same sites. We 
considered that the main contribution comes from the 
matrix element taken between the same eigenstate and 
this is why we considered only the jry) = \a) term 
in the previous equation. One has to note that the 
original integrals in Eq. (j32p are convergent but when 



one breaks them up, then the integral (a| j^ri^— p-|<^) is 
strictly speaking not convergent because one has taken 
one position term y — out of it. The underlying 
assumption which allows us to make this decoupling is 
that, the orbit radius is never allowed to be smaller than 
the effective atomic orbit of the valence state so that the 
cubic singularity does not occur. 

This term gives rise to a contribution in the first 
cumulant. To obtain it, we start from Eq. ([6]) instead 
of Eq. dS]). We substitute Eq. §^ in Eq. §^ and use 
Eq. ([7]) to transform {a\y\/3) to {a\vy\P). The result is 



2^2 



E 




dSa 



1 



d 



^ 4to2c2(47]-££q) 1 1- _ R,^ 



(33) 



(34) 



Again we have used Eq. ([50)1 . Following Datta^^, we will 
assume that, in a crystal, the sum rule of Eq. pop is 



indeed the effective mass. The sum rule (Eq. ((30|) ) will 
be discussed in detail in Section IV Dl Note that if we 
use Eq. (|16p with broadening, then an exact result can 
be obtained by taking the derivative with respect to the 
broadening (^) instead of the energy in Eq. ((33|) . This 
relation will be used to derive Eq. (|36p for the strong 
scattering limit. 



D. The problem of the sum rule of Eq. PP)) 

There is a certain arbitrariness in the use of the sum 
rule of Eq. ((30)) which we should clarify. We note that, 
if we use the standard transformation of Eq. ([7]), then it 
follows, in principle, that when we sum over the entire 
real spectrum of the Hamiltonian, we obtain the free 
electron mass on the LHS. This is simply a consequence 
of the trivial identity jt^x — xpx = —ih with = mvx- 



However, if we follow Datti 



,19 



and use the linear 



response density matrix to compute the acceleration of a 
particle in an electric field, then we obtain the effective 
mass, in the sense of Newton's law, as given by Eq. ()30p . 
Consequently, this gives the absurd result that the 
accelerating particle is always free, irrespective of what 
its initial state is. In effect, the trivial result signifies 
that if we wait long enough, then even a strongly bound 
electron will eventually be free in an electric field. That 
is, this result (mass is free mass), would represent the 
very long time behavior, when the history of the particle 
is irrelevant, and its acceleration in a constant field is 
truly dominated by what happens when it has reached 
its final free state. Dattapi^ concluded that he should 
use a finite band in the evaluation of Eq. pop , and then 
the LHS is indeed the effective mass in the sense of the 
tight-binding band structure, for example. 

The solution of this problem, in general, seems to 
be that, in a transport situation, where electrons are 
injected at one end and absorbed at the other, the sum 
can only run over that part of the spectrum which is 
accessible to the carrier in its lifetime, i.e. for which 
l('^l^2;|/3)P<^(e/3 — £a — fioj) is finite. In weak 
scattering, the particle lives in energy levels near the 
Fermi level which have an effective mass, because it 
relaxes and emits energy to the lattice. In a strong 
scattering situation, the kinetic energy of the carrier 
can be of the same order as the scattering energy 
uncertainty ^. So here, we can relate the sum on the 
RHS of Eq. (1201) directly to the quantum diffusivity D° 
(see Appendix ID)) . We proposes therefore, in the strong 
scattering limit, where Bloch's theorem does not apply, 
to define the effective mass, in the sum rule of Eq. (|30)) 
by the relation 



1 



2Ca 



(35) 
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where Cq is a constant 



^ which carries a 

/3 



sign and is averaged over the band (See Appendix [D] 
for a formal representation). It is of order 1 when the 
energy is near the bottom of the band or in a rapidly 
changing region of the DOS. Wc shall henceforth absorb 
this constant in the definition of an effective diffusivity, 
Da- Equation pS]) is exact f Appendix [Pjl . 

In the random phase limit, Cq is a relatively weak 
function of energy and can be treated as a constant. 
The unit of time is the scattering time. Acceleration 
with strong disorder is therefore drift velocity divided 
by scattering time. With the same definition, in the 
semi-classical limit we therefore have 



1 



d 

dsa V 2m 



(36) 



It follows that a localized initial state |a) has no accelera- 
tion without phonons, i.e. at zero temperature, its effec- 
tive mass in the sense of Eq. ([50)) is infinite. If one evalu- 
ates what dilfusivity one needs to reproduce the electron 
mass using Eq. ([35|) . one has D ~ 1 cm^s~^ which is 
not a small value in a disordered system. We now have 
a way of interpreting terms involving the effective mass 
ma- For example, for strong disorder, expressions of the 
form Ylia ^ ' encountered in Appendix [Bl for 
the second order cumulant can be written as 



D(e)p{e) d 
h de 



1 daje) 
h de 



(37) 



with a{e) = p{£)D(e), where cr(e) is the energy depen- 
dent conductivity with universal scaling properties near 
mobility edges^. 



E. The resistance dependence of the anomalous 
term 

Experimentally one is normally interested in the 
resistance dependence of the AHE and this is the 
criterion used to discriminate between possible models. 
For example, the observed resistance independence of 
the experimental anomalous Hall conductivity in dilute 
magnetic semiconductors is taken as proof that the AHE 
is intrinsic. However, it follows from the present analysis, 
that the apparent scattering-time independence of the 
anomalous Hall conductivity, as derived here in the form 
of Eq. and by Karplus and Luttinger— , and by 

Jungwirth et al7 constitutes a very special limit. The 
scattering time independence of the g-shiit (Eq. 
follows when the matrix elements are dominated by 
interband processes where the energy differences are 
^ ^, (see Ref. [l^ for example). But this implies that 
the basic band structure is Bloch-like, with well-defined 
semiconductor bands. This is obviously not always the 
case, and there can be many situations where disorder 



and band crossings wash out the Kane subband gaps, 
and give rise to arbitrarily small energy denominators 
in which the lifetime directly enters the AHE as well. 
Thus, the g-shift (Eq. ([25]) ) can very well involve the 
conductivity scattering time. 

The assumption of Bloch functions is made by Jungwirth 
et al7 , and the resistance independence of their result 
also rests on the existence of well defined subbands, 
and the sums run over all occupied levels. Thus we 
conclude that for strong spin-orbit scattering and weak 
disorder scattering, the anomalous Hall conductivity 
will not depend strongly on the resistance, and this then 
accounts for some of the experimental observations on 
the AHE^^. In contrast, for weak spin-orbit coupling 
and strong disorder, the anomalous Hall conductivity 
can vary as the scattering time r", < n < 2. The 
Bloch matrix elements^ are not appropriate in the 
strong disorder limit, and the present approach, though 
perturbational, is more appropriate. 



VI. DISCUSSION 

The Kubo formula was expanded in a cumulant expan- 
sion which converges quickly at high temperatures. One 
can also use the series expansion at lower temperatures, 
provided one examines the higher order contributions 
for convergence. Working with the effective mass Hamil- 
tonian, and keeping only the contributions to first order 
in magnetic field allows one to stop the expansion after 
the second order cumulant . The higher order cumulant 
contributions are examined in the Appendices [B] and [Cl 

The second order cumulant derived in Appendix [Bl also 
has a very interesting structure. It gives terms which 
scale as the derivative of the density of states at the 
Fermi level, both for the normal, and the anomalous 
contributions. Previously, Bush and Guentherodt'^*' had 
suggested that the sign of the Hall effect in disordered 
materials scales as the sign of the derivative of the den- 
sity of states at the fermi level (elecron like for increasing 
and vice versa). In the present theory, we have actually 
successfully identified the quantities that determine the 
sign of the Hall effect. Indeed, we have made the very 
interesting observation that the sign of the normal term 
scales, as expressed in the first term, related to the 
derivative of the energy with respect to i?-field. In the 
second cumulant contribution, the sign is related to the 
derivative of the density of states at the Fermi level (at 
low T). Interestingly the two terms have opposite trends 
because an increasing density of states actually gives 
the hole sign, and not as one would intuitively expect 
and as Bush and Guentherodt^S suggest, the electron 
sign. This is a truly remarkable result when applied to 
the strong scattering limit, because normally in weak 
scattering, the first term is dominant and that gives 
the intuitive result from Ref. [sol . But in very strong 
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scattering, for example in the amorphous Hmit, and in 
the region near the mobihty edge, it may happen that 
the second term dominates. If this happens, then we 
have a sign anomaly because the increasing density of 
states at £/ for n-doping gives the hole sign and vice 
versa. This is exactly what is observed in the band edge 
of doped amorphous silicon^. This observation merits a 
more detailed investigation which goes beyond the scope 
of this paper. 

In disordered systems, one can use the CPA (Co- 
herent Potential Approximation)'^^ to describe the 
disordered band structure for exampleSS-^'' and get 
explicit results for the sign of the Hall effect. It turns 
out that there is no simple rule for the sign of the Hall 
coefficient in CPA either, but at the band edges, we 
do indeed have the same behavior as predicted here. 
This has been discussed in detail in Refs. HI and [13. 
In Ref. the case of an impurity band is also considered. 

The disadvantage of the cumulant method is that 
the low temperature limit has to be examined with care 
for convergence. In situations with Bloch symmetry 
where the dominant matrix elements are on the same 
energy shell, this is no problem. But in general, with 
disorder, there is, in the present formulation, in lowest 
order in magnetic field Bz, unfortunately, an infinite 
number of terms. This seems a big problem at first, but 
then it resolves itself. The resolution of the difficulty is 



most obvious when we apply the first order term in the 
limit that the states are localized at £y. We obtain a Hall 
conductivity which is small, but non-zero. This small 
contribution must be canceled by the remaining linear 
terms in the series. Nevertheless the approximation is 
still good because it gives a negligible contribution to 
the Hall current, knowing that the exact result should 
be strictly zero. Eqs. (|B3p and (jCip have terms which 
scale as B^- They are generally smaller than the first 
order contributions we derived in Section |Vl However, 
they involve higher derivatives of the Fermi function and 
can be dangerous to handle at low T. One may infer 
that, if the system has a density of states and scattering 
times which are only weak functions of energy at the 
Fermi level, these higher order cumulant terms are 
negligible. If the density of states is a strong function 
of energy, the expansion will not converge so easily. 
Indeed, near the band edges, there will be mobility edges 
and localized levels which must give rise to a null result 
without phonons, but the null result must be arrived at 
by cancellation of many, albeit, small contributions. 

At high temperatures, the higher order terms cause no 
problem and can be neglected. At any temperature 
we may conclude that a very good approximation is 
obtained by keeping only the first and second cumulants, 
the spin orbit velocity contribution, and the external 
field induced spin orbit term. The final approximate 
formula for the Hall conductivity thus becomes 



f>2 

a 



dSa 



-1 dSa 

'^dB, 



orb 



2m\e 



2in ^ 2ma del [ ^ 



e n \ - 



fl c 
/(ea) 

ma 



E 

a. 



dea ) " dea \ 2ma ) ^ ' 



(38) 



We may call Eq. (|38p the weak-to-intermediate scatter- 
ing Hall conductivity. When we evaluate the cumulants, 
one can use the effective mass Hamiltonian so that the 
periodic potential is not part of the Hamiltonian with 
which the higher order commutators Eq. (|10p arc to 
be evaluated. This means that the only terms which 
contribute above the 2'"^ order cumulant, and which 
scale linearly with Bz^ will be those which depend on 
the disorder and spin-orbit part of the Hamiltonian and 
this makes the approximation of only keeping up to 
second order very accurate. Each of the four additive 
terms of Eq. ((38l) will now be discussed and a simple 
interpretation given. 

The first expression in the bracket of the first term 
has been discussed and is easy to interpret, but it is not 
completely trivial to see that it simply reduces to the 



classical result (Eq. (jlSpl if we used Eq. (fT6|) or Bloch 
functions in the Kubo formula with a constant lifetime. 
In the pure quantum limit it gives the Streda result but 
it is actually the normal Hall effect. The way this term 
should be handled depends on the problem in question. 
In weak scattering it again gives Eq. ([T5)) . with a free 
electron mass. In the tight-binding representation, one 
can evaluate it using second order perturbation theory 
in the magnetic field dependent term in the Hamiltonian. 

The second part of the first bracket involves the g- 
shift of the delocalized levels above the mobility edge. 
In comparison, the localized g-shift is negligibly small. 
In the weak scattering limit, the g-shift can be evaluated 
using the Kane-Luttinger wavefunctionsiS. This is 
given by Roth et aP^. In an incoherent situation, a 
disordered system with no Bloch symmetry and strong 
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scattering , we should use rric 



as the effective mass. spin-orbit magnetic field. 



For the second term, the one involving a second 
derivative of the Fermi function, one can, using the 
nia — -jj- approximation, obtain integral products of 
the type a{e) = p{e)D{e) as shown in Eq. ([57]) which are 
energy dependent conductivities, and which obey well 
known universal scaling relations near the mobility edges. 

The third term was shown to be due to the effect 
of the external-field-induced spin-orbit energy (Sec- 
tion |VB]). With TTia — m*, this term reduces to the 
same form as the so called side jump contribution. 
Our theory shows that it can be enhanced via a small 
effective mass and even extended to apply to strong 
scattering via Eq. ([55)1 . 

Considering the second and third terms, we pro- 
pose that, to a good approximation, in most situations 
where Bloch functions cannot be used, we may replace 
nia = -jj-, once the sign has been determined via 
Eq. (l30|) . This accounts for localized states if any are 
present, because Da is zero. 



The last term in Eq. ([55]) is due to the spin-orbit 
contribution to the velocity operators via the internal 
potentials from Eqs. ^ and ([5]). This also vanishes 
for localized states. Above the mobility edge, this 
contribution is comparable to the second expression 
in the bracket of the second term. This can be seen 
using the strong scattering correspondence (Eq. (l36p V 
The last term in Eq. (|38p then reduces to a term which 
resembles the usual skew scattering terni^^ provided we 
interpret the skew scattering rate for Coulomb spherical 
potentials as 



1 



he 



r-R„| 



a 



(39) 



Thus, for disordered systems, the last term of Eq. ([38 
term can be written 



_ske'w 



dSa 



T mDr, 



(40) 



Note that the sum runs over all the potentials in the 
lattice. So the skew scattering and g-shiit terms include 
both extrinsic and intrinsic contributions. Note that 
very recently, Chudnovsky-S^ using a different approach, 
and for spin Hall effect, also obtained a term where all 
the potentials are included. In the form of Eq. (|40)) . 
appropriate for disordered systems with no Bloch 
symmetry, the Bloch enhancement does not appear. If 
we neglect the host spin-orbit coupling and only include 
the impurities, then both f;-shift and skew term are, 
by definition, extrinsic contributions, and the mass 
in Eq. ([M)) is necessarily the effective mass. It is the 
effective mass particle which generates the impurity 



Let us now re-examine the question of the theoret- 
ical resistance dependence of the AHE. From the above 
analysis we note that this all has to do with the way 
we treat the matrix elements, and at what stage we 
introduce incoherence and lifetime. This can already be 
seen in the first term, which can be treated as a quantum 
effect, as in Stredapi^ or in the semi-classical limit. The 
same is true for the spin-orbit terms. Thus, if we keep 
to the notion of effective mass, we have the quantum 
result. If we use the transformation of Eq. ([55]) . which 
involves the diffusivity, then we have the connection 
with conductivity. As an example, we can take the first 
two, dominant, terms in Eq. (|38p and write them, using 
the definition of the g-shift given by Roth et al?^ (we 
assume A and Eg to still be defined), in the incoherent 
limit as: 



eB,T 



-T. 



dfjeg) 

dSa 



m 



Da , 

m— 1 , „ 

h / 3E, 



A 



2A 



(41) 



where we have used Eq. p5p in the g-shift term. 



Now, we see that what was a pure quantum term with 
effective mass, has, in this limit, become a term which 
depends on the diffusivity. The AHE Hall coefficient can, 
it seems, change from a linear scaling with resistance 
in Ea. (|^T|) to a resistance squared (relaxation time 
squared) behavior if becomes m~^, as is observed in 
diluted magnetic semiconductors^^. The experimentally 
observed Ra oc p^^ dependence implies, therefore, that 
the effective mass concept remains valid in these systems. 



We have focused our attention on the anomalous 
contributions, and how they undergo a transformation, 
in going from the weak scattering to the strong scattering 
limit. A similar change must occur for the first term of 
Eq. (PT|) . Here too, we must replace the Drude a^x by 



de - 



dfje) 
de 



p{e)D{e) 



and the band mobility term 



by 



erB, D 
— = eB^ — , 

m* (e) 



(42) 



(43) 



at low temperatures. The notation in terms of diffusivity 
D is valid in the weak scattering limit too, but now one 
can see what happens as we approach the mobility edge 
using standard localization theory and a{e) — p{e)D{e). 

We have shown that the AHE Hall coefficient 
can vary with resistance p"^ with 1 < n < 2 depending 
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on the degree of coherence. The essential point seems to 
be the way one treats the sum rule given by Eq. ([50)) . 
This was already a problem for Dattaii^ in his analysis 
of the Hall effect in his 1980 paper. Eq. ([50)1 may be 
treated as the well known f-sum rule, but it is clearly not 
realistic to sum over an infinite excited state spectrum, 
without taking into account the finite lifetime of the 
states. The question then becomes: when is the sum 
on the RHS of Eq. ([30| a-nd when is it closer to 
which is the sum evaluated in the semi-classical 
limit, and also derivable from (Q;|a;^|Q;) ~ DaT? The 
interesting and important side of this last result, is that 
it gives the correct null result for transport in localized 
states, and therefore must be the correct approximation 
near the mobility edges. 

We have derived a side-jump-like contribution to 
the Hall current using linear response to the (Rashba) 
spin-orbit interaction^^ which is produced by motion 
in the external potential. Our result is closely related 
to the Rashba current— and is in principle, exact, up 
to the evaluation of the sum rule (Eq. ([30)) ). The sum 
rule has to be treated with care in both the coherent 
and the strong scattering limit. In the incoherent limit, 
and specially near mobility edges, we have the relation 
nia = which then gives a very appealing result, 
namely 



'xy 



(44) 



This gives a vanishing side-jump-like Hall effect (Rashba) 
contribution from those states in the localized part of 
the band. 

Adding Eqs. ([44| and Eq. ([4T)) we have the strong 
scattering Hall conductivity limit. As in the QHE, 
every conducting state contributes. However, the full 
enhancement of the spin-orbit coupling (see Ref. [13 and 
Ref. [l^ for example) does not appear in our expres- 
sion. The most interesting aspect of the side-jump-like 
(Rashba) term (Eq. ([11])), is that if the Fermi level 
is in a region of localized states, say near the top of 
the band, then this term dominates the Hall current 
since all other terms are Fermi level terms, and vanish 
at low temperatures. It would give us (using the free 
mass) Jy ^ lO~'^^Ndei{(Jz)Fx A/m^ (Ndei = delocalized 
electrons density) which is, of course, a small but rigid 
current, analogous to the Quantum Hall current. 



VII. CONCLUSION 

The principal result of the cumulant method developed 
in this paper is given by Eq. ([38)) . It is an expression for 
the Hall current which allows for disorder and spin-orbit 
coupling. For a typical doped semiconductor, we might 
expect the following contributions to the Hall effect: 



1. The normal process (first term of Eq. ([IT))). 

2. The intrinsic AHE side jump like process (Rashba 
term, Eq. 



3. The intrinsic AHE Karplus-Luttinger process (sec- 
ond term of Eq. ([iT)) ). 



4. The skew scattering process due to impurities 
(Eq. ([10])). Here, the basic spin-orbit coupling has 
the effective mass and not the free mass. 

The skew term (Eq. (HUl) ) is only relevant if both the 
enhancement and the basic spin-orbit coupling are large 
and the concentration of impurities high enough (see 
Ref. [1^ for estimates). The side-jump term is only impor- 
tant if the effective mass is small or the Fermi level is in 
a region of localized states. In most situations with weak 
disorder, we expect the Karplus-Luttinger intrinsic term 
to dominate the AHE. In the limit of strong scattering, 
when the use of Bloch functions and uniform lifetimes 
is no longer valid, we may replace the effective mass in 
Eq. ([38]) using the concept of quantum diffusivity i.e. put 
TO~^ — > and use derivatives, as explained in the text. 

For magneto-optical, Faraday angle measurements, 
for example, we need the frequency dependent Hall 
conductivity. One can show that in weak scattering, 
where we use the effective mass for Eq. (I5(I1) . we add 
a factor i^^lj^yi to the four main Hall conductivity 
results listed above. In strong scattering, where we use 
the diffusivity to describe Eq. ([5(I)) . D{e) is replaced by 
D{e,Lj). 
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APPENDIX A: DERIVATIVE OF THE 
HAMILTONIAN WITH RESPECT TO THE 
MAGNETIC FIELD 

With Ai = {Q,BzXi) and Py = mVy, the kinetic part 
of the Hamiltonian is 



0\2 



2m 



(Al) 



Thus, with H ^ T + V + Hso ~ ^B, J2^ <^ we can 
write, remembering that without the spin-orbit term in 
the velocity operator v't 



0,i _ eB^ 

^ '^y m -^J' 



dB, 



eXiV 



V 



(A2) 
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Thus, 



dH 



-e{a\xVy 



(A3) 

We now want to transform the LHS of Eq. (jA3p to ob- 
tain Eq. (I13|) . To do so, we just apply the weU known 



dH 



Hellman-Feynman theorem and we obtain i^a 

and thus Eq. Remember that Vy in the equa- 

tion does not contains the spin-orbit term of Eq. ([5]). 



APPENDIX B: THE SECOND ORDER TERMS 
IN THE CUMULANT EXPANSION 

Consider now the second order contribution to the Hall 
conductivity. This can be written 



.{2} 



— Y^{a\x\P){mvy-VvH\a) 

a, 13 



del 



(Bl) 



where, as previously, Vy = — dropping the 

spin-orbit current operator. After evaluating the com- 
mutator we are left with three terms 

_eBz n ihdV ih 



ih- 



'z 
-?; 



m dy 



-Px 



^nCTz 



(B2) 



Only the first and third term give significant contribu- 
tions in this order. The first term is 



2in 



^2 h^eB 



del 



2m 



m 



2rar 



del ■ 



(B3) 



where we have used the sum rule of Eq. ([50)1 restricted 
to the active energy band. For Bloch states restricted to 
the 8-Kane bands with spin-orbit coupling, the RHS of 
Eq. ((30|) reproduces the 8-bands k • p calculated effective 
masses. To conclude, we note that the cumulant expan- 
sion is a powerful method when the bandwidth is very 
narrow and < fc^T. At low temperatures, when we use 
Vila = TTL* , we obtain for Eq. (|B3p 



/2} 
'xyA 



h f heBz \ dp[e) 



2\ m \ 2m* 



de 



(B4) 



Note that the cumulant expansion is not purely an ex- 
pansion in powers of B^. There is a term linear in B^ in 
almost every order. Using the effective mass method, the 
linear term in Bz becomes an expansion in powers of the 
Vdisorder / e f . The Convergence is problematic near band 
edges where we have localization. We neglect the term 
given by the second term of Eq. (|B2p . The spin-orbit 
term, given by he third term of Eq. (|B2[) . becomes 



xy,3 



e2ft2 



1 



2in ^ 2ma 



d'fjeq) 

' del. 



(B5) 



Again, at low T (neglecting other energy dependence un- 
der the integral), we can use 



Ct 



del 



dp 

de 



(B6) 



again, we have the 



If we assume that ^1 
same structure as before in Eq. psp . this time with 
^5 = ■^^T^- With the numerator of order 10"'* eV, this 

term corresponds to an effective g-shift of 10"^ which 
is therefore smaller than the first order cumulant of its 
type. But, in general we remind the reader that the cu- 
mulant expansion, as is also true for the configurationally 
decoupled Kubo formula, is not accurate near the band 
edges for reasons of Anderson localisation. When used in 
this region it can only apply above the mobility edge. 



APPENDIX C: THIRD ORDER CUMULANT 

In the third order cumulant, if we neglect terms of 
order B^ and X^, we have only the terms 



xy 



—y 



1 d^fiea) 





The terms linear in Bz and spin-orbit coupling always 
go in pairs, the spin-orbit term acting like an effective 
magnetic field. The higher order terms form an infinite 
series with products involving derivatives of the lattice 
potential. These terms are small at high temperature 
and renormalize the first order linear terms in Bz and 
A„. In the effective mass approach, the lattice potential 
no longer appears in the Hamiltonian so ^(r) appearing 
in Eq. (jCl|) is due to impurities or disorder. 



APPENDIX D: EFFECTIVE MASS AND 
DIFFUSIVITY IN THE STRONG SCATTERING 
LIMIT 

Since the excited states can decay and must have a 
finite lifetime, we can write the sum of Eq. ([50)1 as 

1 ,^ J(aK|/3)|2(e^-e,) 
— -— = 2 / depp{efi) — -y. (Dl) 

""^^"^ (e/3-£.)2+(A) 

The quantum diffusivity is defined as 

Da = 'hY,\Hvxm^Ke^~ea). (D2) 
/3 
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Then it follows that we can write 
1 



2 / dsfil p{ep) 

2CaD{£a) 



A 



(D3) 



where is a constant, defined by Eq. (jD3p . Its value and 



sign depends on the density of states variation. One can 
also see this by doing an integration by parts using the 
product of the two functions in brackets. The constant 
carries a sign (as does the effective mass) which depends 
on the energy. We can rewrite it as Cq, — {{ep — Eq)-^)^, 
where the average is to be taken with the diffusiv- 
ity weighting function under the integral as defined by 
Eq. (|D3p . This equation is still exact. In the random 
phase approximation, kinetic energies are of the same 
order as the energy uncertainty and Cq, is —1 or +1 
depending on whether we are near the top or the bottom 
of the band. 
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